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ORIGINAL ARTICLE 

Effects of genetic and early environmental risk factors for 
depression on serotonin transporter expression and 
methylation profiles 

M Wankerl 1 , R Miller 1 , C Kirschbaum 1 , J Hennig 2 , T Stalder 1 and N Alexander 1 

The serotonin transporter (SERT) gene-linked polymorphic region (5-HTTLPR) has been implicated in moderating the link between 
life stress and depression. However, respective molecular pathways of gene-environment (GxE) interaction are largely unknown. 
Sustained alterations in SERT gene expression profiles, possibly mediated by epigenetic modifications, are a frequent correlate of 
depression and may thus constitute a putative mediator of GxE interaction. Here, we aimed to investigate joint effects of 5-HTTLPR 
and self-reported environmental adversity throughout the lifespan (prenatal, early and recent stress/trauma) on in vivo SERT mRNA 
expression in peripheral blood cells. To further evaluate whether environmentally induced changes in SERT expression are 
mediated by epigenetic modifications, we analyzed 83 CpG sites within a 799-bp promoter-associated CpG island of the SERT gene 
using the highly sensitive method of bisulfite pyrosequencing. Participants were 133 healthy young adults. Our findings show that 
both the 5-HTTLPR S allele and maternal prenatal stress/child maltreatment are associated with reduced in vivo SERT mRNA 
expression in an additive manner. Remarkably, individuals carrying both the genetic and the environmental risk factors exhibited 
32.8% (prenatal stress) and 56.3% (child maltreatment) lower SERT mRNA levels compared with those without any risk factor. Our 
data further indicated that changes in SERT mRNA levels were unlikely to be mediated by DNA methylation profiles within the SERT 
CpG island. It is thus conceivable that the persistent changes in SERT expression may in turn relate to altered serotonergic 
functioning and possibly convey differential disease vulnerability associated with 5-HTTLPR and early adversity. 
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INTRODUCTION 

Research has consistently implicated the joint contribution of 
genetic and environmental risk factors in the pathogenesis of 
major depression. 1 The most prominent example refers to the 
debate on whether a 43-bp insertion/deletion polymorphism (5- 
HTTLPR) in the serotonin transporter gene 2 moderates the 
association of life stress and depression. 3-5 Whereas numerous 
studies observed increased disease vulnerability in carriers of the 
5-HTTLPR short (S) variant upon exposure to environmental 
adversity, 4,6 recent meta-analyses have triggered an active 
controversy about whether this finding holds up. 7,8 Therefore, 
exploring systemic and molecular mechanisms underlying gene- 
environment (GxE) interactions seems to be of major importance 
to further advance this debate. 1,9 

On a systemic level, experimental studies investigating biolo- 
gical quantitative traits strongly support the hypothesis of 
elevated stress sensitivity in S allele carriers. 4 Among other 
biological alterations, the S allele has repeatedly been associated 
with elevated amygdala activity 10,11 and increased Cortisol 
secretion in response to a variety of aversive/stressful 
stimuli. 12,13 However, little is known about the molecular path- 
ways mediating disease vulnerability. One hypothesis is that GxE 
interaction already takes place at the very early level of gene 
expression, in a way that 5-HTTLPR and life stress jointly convey 
stable changes in serotonin transporter (SERT) expression. In line 



with this, altered SERT expression profiles may constitute a 
putative mediator of GxE interaction as they have been commonly 
observed in depressed patients 14-16 and stress-sensitive Rhesus 
macaques} 7 '™ 

The functional effects of 5-HTTLPR have been widely docu- 
mented by in vitro studies, indicating that the S allele is associated 
with reduced SERT gene (SLC6A4) transcription in lymphoblast cell 
lines 2,19 and decreased serotonin uptake in platelets. 20,21 In 
addition, transcriptional efficiency of the SERT gene was found 
to be influenced by an A/G single-nucleotide polymorphism 
(rs25531) located upstream of the 5-HTTLPR promoter variant 
within the greater repeat structure. 19 This has led to the 
distinction between the variants S, L A and L G , (tri-allelic classifica- 
tion of 5-HTTLPR) with the latter one being functionally similar to 
the S allele. 19 In contrast to in vitro studies, results obtained in vivo 
appear to be less conclusive regarding allele-specific SERT mRNA 
expression in peripheral cells 22,23 and SERT availability in the 
human brain. 24-26 Environmentally-induced changes in SERT 
expression may partly account for the observed 
inconsistencies 27 but have been sparsely addressed in humans. 
Animal studies in various species provide first evidence that 
exposure to early adversity correlates with decreased SERT mRNA 
levels in the brain 28,29 (but also see Gardner et al. 30 ) and in 
peripheral cells. 31 These long-term changes in SERT expression 
patterns may result from stable epigenetic modifications such as 
DNA methylation. 32 Recent studies using peripheral blood cells 
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found that increased methylation levels within a 799-bp 
promoter-associated CpG island in 5LC6A4 associate with both 
lower SERT mRNA levels 23,33-35 and exposure to childhood 
trauma, 23,36-39 in some studies dependent on 5-HTTLPR 
genotype. 40 Such peripheral measures of gene expression and 
DNA methylation profiles have been increasingly recognized as 
informative biomarkers in psychiatric research. 41-43 Most impor- 
tant for the present study, epigenetic and transcriptional changes 
in response to environmental adversity appear to be system- 
wide 44,45 and can thus potentially be tracked in easily obtainable 
blood cells. 

The present study aimed to investigate joint effects of 5-HTTLPR 
and environmental adversity across different developmental 
stages (prenatal, early, recent trauma/stressors) on peripheral 
SERT mRNA expression and DNA methylation within the 
promoter-associated SERT CpG island. As a potential molecular 
pathway of GxE-mediated disease vulnerability, we expected to 
find lowest SERT mRNA and highest SERT methylation levels in 
individuals carrying both the genetic (S allele) and environmental 
(stress/trauma) risk factor. Studies investigating long-term tran- 
scriptional signatures of early adversity implicitly assume that 
gene expression profiles are characterized by a trait-like compo- 
nent with substantial differences between individuals. Prior 
studies have generally confirmed a considerable intraindividual 
stability of genome-wide gene expression patterns over hours and 
months; however, they have also highlighted that stability varies 
across individual transcripts. 46,47 As intra- and interindividual 
variation in human SERT mRNA expression has, to the best of our 
knowledge, not explicitly been evaluated yet, we further 
conducted a small pilot study investigating respective patterns. 

MATERIALS AND METHODS 

Pilot study 

For the assessment of intra- and interindividual variation in SERT mRNA 
levels, we obtained SERT mRNA expression day profiles in eight healthy 
individuals (four females, mean age: 25.0 ±3.3 years) at two test days 
separated by 1 week. On each day, seven blood samples were drawn from 
an indwelling cannula (every 2 h from 0800-2000 hours) into PAXgene 
blood RNA Tubes (PreAnalytiX, Qiagen, Hilden, Germany). As we were 
interested in normal diurnal fluctuations of SERT mRNA expression, no 
restrictions were imposed on participants regarding food intake, work 
flow etc. 

Main study: sample and procedure 

We recruited participants aged 18-30 years via newspaper advertisement 
and flyers. Only healthy Caucasian participants who were native German 
speakers were included in the study. After a structured telephone 
interview that served as a first screening for exclusion criteria (for example, 
major health issues), 155 individuals were invited for the main screening 
and testing session. During this session, the Diagnostic Interview for 
Psychiatric Disorders — short version (Mini-DIPS 48 ), a structured interview 
assessing point and lifetime prevalence of axis I disorders based on DSM IV 
criteria, was conducted. Furthermore, participants completed a compre- 
hensive checklist on chronic physical diseases (for example, cancer, 
diabetes, heart diseases, asthma and epilepsy) and medication intake (for 
example, psychotropic drugs). Any current or past mental and/or physical 
disease as well as medication intake and pregnancy were defined as 
exclusion criteria. Participants who passed this screening procedure (final 
sample size: A/=133, 63 females) were asked to fill in a set of 
questionnaires on early and recent life stress/trauma. Data on prenatal 
stress were obtained within a subsample of 85 participants of whom their 
mothers agreed to fill in a questionnaire on maternal stress/trauma during 
pregnancy. At the end of the session, blood samples were drawn into EDTA 
tubes (Sarstedt, Numbrecht, Germany) for DNA and PAXgene blood RNA 
Tubes (PreAnalytiX, Qiagen) for RNA extraction and stored at - 20 °C for no 
more than 6 months. The pilot and the main study were conducted in 
accordance with the Declaration of Helsinki and were approved by the 
ethics committee of the Technische Universitat Dresden. Participants 
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provided written informed consent and received a monetary reward for 
participation. 

Assessment of prenatal, early and recent life stress/trauma 
In order to assess prenatal stress/trauma, mothers of participants 
completed the NeuroPattern-PreVpostnatal-Stress-Questionnaire, which 
retrospectively records pre-, peri- and postnatal adverse events. 49,50 The 
NPQ-PSQ is part of a translational diagnostic tool (NeuroPattern) for stress- 
related disorders 49 and assesses maternal stressful/traumatic events during 
pregnancy, such as death or life-threatening illness of a close relative, 
divorce, lack of social support, relationship conflicts, high workload and 
financial constraints in a yes/no format. Participants of mothers reporting 
at least one stressful/traumatic life event during pregnancy were assigned 
to the prenatal stress group. We further applied the short Form of the 
Childhood Trauma Questionnaire (CTQ), 51 ' 52 a widely used retrospective 
measure of child maltreatment with high internal consistency, reliability 
and criterion validity. 51 Participants were classified as traumatized when 
CTQ scores exceeded a moderate to severe cutoff score 51 in at least one of 
the five CTQ trauma categories (emotional abuse: >13, physical abuse: 
>10, sexual abuse: >8, emotional neglect: >15 and physical neglect: 
>10). In addition, recent stress exposure was assessed using the Life 
Stressor Checklist— Revised (LSC-R 53 ' 54 ). The LSC-R is a 30-item self-report 
measure with good psychometric properties 54 assessing traumatic and 
stressful life events (for example, physical/sexual assault, death of a 
relative, serious accidents/diseases, abortion) in a yes/no format. Partici- 
pants reporting at least one stressful/traumatic life event within the past 5 
years were assigned to the recent stress/trauma group. 

5-HTTLPR genotyping 

DNA was extracted from EDTA whole blood using a standard commercial 
extraction kit (High Pure PCR Template Preparation Kit; Roche, Mannheim, 
Germany) in a MagNA Pure LC System (Roche). Participants were 
genotyped for 5-HTTLPR/rs25531 according to a previously published 
protocol 55 

Quantitative real-time PCR 

Real-time quantitative PCR was performed by Varionostic GmbH (Ulm, 
Germany; http://www.varionostic.de) on the LightCycler 480 I (Roche) 
using the SensiFast Sybr green mix from Bioline (Luckenwald, Germany). A 
detailed protocol and primer sequences are provided in Supplementary 
Information 1. Glyceraldehyde-3-phosphate dehydrogenase and b-actin 
were applied as references for expression. After the mean calculation and 
delta Cq generation (reference mean-target gene value) values were 
plotted according to the 2~ 2AACt method 56 to form relative expression 
level on a linear scale. Values reflect fold changes in gene expression 
normalized to both endogenous reference genes. All RT-PCR analyses were 
performed in duplicates. 

Bisulfite pyrosequencing 

We analyzed quantitative methylation of 83 CpG sites within a 799-bp 
promoter-associated CpG island in SLC6A4 (GenBank accession number: 
NG_011747). Methylation analysis by bisulfite treatment of genomic DNA 
from EDTA whole blood samples and subsequent pyrosequencing was 
performed by Varionostic GmbH. Sequencing was performed on the Q24/ 
ID System and percent methylation at each CpG site was quantified using 
the PyroMark Q24 software (Qiagen). A detailed protocol with amplicon 
and sequencing primers is provided in Supplementary Information 2. The 
percentages of methylation values, which passed quality control, were 
>95% for each individual CpG site. 

Statistical analyses 

Analyses were conducted using SPSS (Version 21.0, IBM, Chicago, IL, USA) 
and R (R Core Team, 2013). Within the pilot study, intraclass-correlation 
coefficients (ICCs) for the factors 'time' and 'subject' were calculated to 
assess intra- and interindividual variation in SERT expression according to 
the following model: mRNAj = j8 SU bject+ftime+£i/ witn £ i ~ A/(0,of es )- 57 Large 
values indicate that a major fraction of variance in SERT expression can be 
explained by the respective factor. Furthermore, SERT mRNA area under 
the curve with respect to ground (AUC G ) values were calculated as an 
integrated measure of total SERT mRNA output according to the 
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trapezoidal formula. 58 Pearson correlations of SERT mRNA AUC G values for 
the two test days were further calculated. 

In the main study, X 2 tests for dichotomous and analyses of variance for 
continuous measures were used to examine group differences regarding 
demographic characteristics. Effects of 5-HTTLPR and stress/trauma-related 
measures on SERT expression, mean as well as principal component 
methylation across the CpG island and methylation at individual CpG sites 
(Bonferroni-corrected), were tested by general linear models. We have 
further specified a general linear model incorporating a GxE interaction 
term in addition to main effects of genotype and stress/trauma. This model 
was compared with a baseline model assuming additive effects only to 
evaluate whether GxE interaction explains incremental variance. Regarding 
genotype as between-subject factor, analyses were conducted both 
according to the bi-allelic (SS versus SL versus LL) and the tri-allelic 
classification by comparing the load of high (L A ) and low (S, L G ) expressing 
alleles (SS,SL G ,L G L G versus SL A ,L G L A versus L A L A ). Furthermore, mediation 
analyses were conducted assessing the relative contribution of SERT 
methylation as a mediator of the association between environmental 
adversity and SERT mRNA expression. 59 Pearson correlations were 
calculated to further evaluate associations between stress/trauma scores 
and the dependent variables. To identify potential confounds, we tested 
for associations of sex, oral contraceptives use, body mass index, age and 
smoking status with our dependent variables as these factors have 
previously been found to influence gene expression and methylation 
patterns (for example, the well-documented effects of sex hormones on 
gene expression and DNA methylation 60,61 ). These analyses were 
conducted by means of independent t-tests for dichotomous and Pearson 
correlations for continuous variables. 

RESULTS 

Pilot study: intra- and interindividual variation in SERT mRNA 
expression 

Figure 1 displays SERT expression levels at the 14 different 
measuring times for each participant and the mean SERT 
expression patterns across the 2 test days. Our findings revealed 
that only a small portion of variance in SERT mRNA expression 
patterns is bound by the factor time (ICC time = 0.05), whereas ICCs 
indicate a moderate to high between-subject variance (ICC subject = 
0.60). We further observed a high intraindividual day-to-day 
stability of SERT expression patterns as indicated by a strong 
correlation (r = 0.89, P = 0.003) between SERT mRNA AUC G across 
the 2 test days. Our pilot study thus indicates a substantial trait 
component regarding SERT mRNA expression, given that a higher 
portion of variance in SERT mRNA expression can be explained by 
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the factor subject (60%) than by other factors (ICC time +residual 
variance = 40%). 

Main study: sample characteristics 

Demographic and stress/trauma-related sample characteristics are 
depicted in Table 1. There was no significant deviation from 
Hardy-Weinberg equilibrium using bi-allelic (xf 1) = 0.72, P = 0.40) or 
tri-allelic (xf 3) = 2.95, P = 0.40) classification of 5-HTTLPR. Genotype 
groups did not differ regarding the number of reported prenatal, 
early and recent stressors/trauma (all P-values ^ 0.08, Table 1) or 
with respect to sex, age, body mass index, smoking status and oral 
contraceptive use (all P-values^O.21, Table 1). When participants 
were assigned to the prenatal, early and recent stress/trauma 
groups, no significant differences between the respective 'stress/ 
trauma' and 'no stress/trauma' groups regarding any of the these 
variables were found (all P-values^0.25), except of a higher 
number of smokers in the 'early trauma' compared with the 'no 
early trauma' group (64.7% versus 30.2%, xf l} = 7.81, P = 0.01). 

Effects of 5-HTTLPR and stress/trauma-related variables on SERT 
mRNA expression 

Table 2 presents effects of 5-HTTLPR and stress/trauma-related 
variables across the lifespan on SERT mRNA expression. SERT 
mRNA levels were found to be unrelated to sex, age, body mass 
index, smoking status and intake of oral contraceptives (all P- 
values^0.16). As expected, we observed significantly lower SERT 
mRNA levels in S allele carriers compared with individuals carrying 
two copies of the L allele (F 1/131 =5.14; P = 0.03; q 2 = 0.04). Similar 
results were obtained with tri-allelic classification of 5-HTTLPR 
(F 1J31 =6.67; P = 0.01; a? 2 = 0.05). On a genotype level (SS/SL/LL), 
analyses of variance revealed a borderline significant effect of 5- 
HTTLPR on SERT mRNA levels (F 1J30 = 2.71; P = 0.07; A] 2 = 0.04), 
which reached significance with tri-allelic classification 
(F 1/130 = 3.31; P = 0.04; A] 2 = 0.05). 

Regarding stress/trauma-related variables, we observed a 
significant effect of prenatal and early stress/trauma on SERT 
mRNA levels. Individuals of mothers reporting at least one major 
stressful/traumatic life event during pregnancy were found to 
have lower SERT mRNA levels compared with those without 
maternal prenatal stress (F 1;83 = 5.54; P = 0.02; q 2 = 0.06). Further- 
more, we observed a negative correlation between the number of 
prenatal maternal life stressors/trauma and SERT mRNA levels 
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Figure 1. (a) SERT mRNA expression levels (mean ±s.e.m.) calculated with the 2 _2AACt method at the 14 different measuring times displayed 
for each participant separately. Intraclass-correlation coefficients (ICCs) for the factors 'time' and 'subject' were calculated to assess intra- and 
interindividual variation in SERT mRNA expression patterns. ICCs can vary between 0 and 1. Large values indicate that a major fraction of 
variance in SERT expression can be explained by the respective factor, (b) SERT mRNA expression levels (mean ±s.e.m.) calculated with the 
2 -2AAct method f or the 2 test days. 
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Table 1. Demographic and trauma-related sample characteristics (mean and s.d.) divided by 5-HTTLPR genotype 



5-HTTLPR Genotype 





Total (n = 133) 


S allele carriers (r\=86) 


L homozygotes (n=47) 


P 


Sex (% females) 


47.4% 


43.0% 


55.3% 


0.21 


Age (years) 


23.82 ±3.05 


23.90 ±3.1 4 


23.68 ±2.89 


0.70 


Body mass index 


22.47 ±2.20 


22.47 ±2.1 2 


22.49 ±2.37 


0.96 


Smokers (% yes) 


35.3% 


33.7% 


38.3% 


0.71 


Oral contraceptives (% yes) 


25.6% 


22.1% 


31.9% 


0.22 


Prenatal trauma/stress (% yes, NPQ-PSQ)* 


43.5% 


42.6% 


45.2% 


0.83 


Early trauma/stress (CTQ Score) 


33.02 ± 7.92 


33.02 ±8.1 3 


33.00 ±7.61 


0.99 


Recent trauma/stress (LSC-R Score) 


1 .43 ± 1 .35 


1.45 ±1.32 


1.38 ±1.42 


0.08 



Abbreviations: CTQ, Childhood Trauma Questionnaire; 51 LSC-R, Life Stressor Checklist-Revised; 53 NPQ-PSQ, PreVPostnatal Stress Questionnaire. 49 *n = 85. 



(r = -0.27, P = 0.01). A similar reduction in SERT mRNA expression 
was seen in individuals reporting early traumatization 
(Fi, i3i =4.15; P = 0.04; q 2 = 0.03) when comparing the groups of 
participants with (mean CTQ score: 48.88 ±4.38) and without 
(mean CTQ score: 30.69 ±8.54) a history of child maltreatment. 
Within the overall sample, CTQ scores were inversely correlated 
with SERT mRNA levels (r = -0.18; P = 0.04). In contrast, recently 
experienced stress had no effect on SERT mRNA levels 
(Fi,i3i = 1.17; P = 0.19). Our findings further indicated that the S 
allele and prenatal/early adversity associate with decreased SERT 
expression in an additive manner as lowest levels were found in 
individuals carrying both the genetic and the environmental risk 
factor (prenatal stress: F 2 , 82 = 3.14; P = 0.05; q 2 = 0.07; early trauma: 
F 2 ,i3o = 4.55; P = 0.01; q 2 — 0.07; Figure 2). Specifically, those 
individuals were found to have 32.8% (prenatal stress) and 
56.3% (child maltreatment) lower SERT mRNA levels compared 
with the group without any of the two risk factors. Including a GxE 
interaction term in the general linear model did not incrementally 
increase the portion of variance explained by an additive model 
(prenatal stress: F 181 =0.86; P = 0.36; early trauma: F ll2 9<0.01, 
P = 0.99). 

Effects of stress/trauma-related variables on SERT DNA 
methylation profiles 

We next investigated whether effects of prenatal and early life 
stress/trauma on SERT mRNA expression levels are mediated by 
DNA methylation profiles within a 799-bp promoter-associated 
CpG island in SLC6A4 (Figure 3a). The mean methylation values 
averaged across the entire CpG island did not differ as a function 
of sex, age, body mass index and smoking status (all P-values^ 
0.14). The use of oral contraceptives was associated with 
significantly increased mean SERT methylation (F 161 =4.20; 
P = 0.05; q 2 — 0.06) and was thus included as a covariate in 
subsequent analyses. 

Our results revealed no significant effect of prenatal (F 182 = 0.1 7; 
P = 0.68), early (F 1J30 = 0.03; P = 0.86) or recent (F 1;130 = 0.73; 
P = 0.39) life stress/trauma on the mean SERT methylation levels 
(Table 2). Furthermore, the mean SERT methylation levels did not 
differ as a function of 5-HTTLPR (all P^0.50, Table 2). Regarding 
functional relevance of SERT methylation profiles, we observed no 
significant correlation between the mean methylation levels and 
SERT mRNA expression (r = 0.10, P = 0.23). The latter results 
indicate that effects of prenatal/early adversity on SERT mRNA 
expression are unlikely to be mediated by overall SERT methyla- 
tion. For the purpose of completeness, we have additionally 
calculated mediation analyses, which overall confirmed this 
presumption (indirect effects: prenatal stress: R = - 0.00003 [CI: 
-0.00340, 0.00277], early trauma: 5 = -0.0001 4 [CI: -0.00468, 
0.00319]). 



As absolute levels and interindividual variation in methylation 
substantially vary across the SERT CpG island (Figure 3b), we 
further conducted exploratory analyses on the level of individual 
CpG sites. First, we screened the entire CpG island for sites related 
to decreased SERT mRNA levels. Methylation levels at 10 out of the 
83 CpG sites investigated were associated with lower SERT mRNA 
expression (all P-values < 0.05 uncorrected, Figure 3a); however, 
only for CpG9 this association remained significant after correcting 
for multiple testing (r = -0.34, P< 0.001). We further observed a 
considerable overlap of site-specific methylation associated with 
SERT mRNA expression and prenatal stress/trauma. Individuals 
exposed to maternal prenatal stress were found to have higher 
methylation levels at CpG2, CpG9, CpG29 and CpG30 (all P-values 
< 0.05) compared with those without. The latter effect remained 
significant for CpG30 (F lf83 = 1 1.81, P = 0.001, q 2 = 0A6) and 
by trend also for CpG9 (F 1#83 = 8.24 # P = 0.005, a? 2 = 0.09) 
after Bonferroni-adjustment (Figure 3a). We further observed 
no associations between early and recent life stress/trauma or 
5-HTTLPR and site-specific SERT methylation levels (all 
P-values ^0.39). 

For completeness, we have further investigated associations of 
prenatal, early and recent life stress/trauma with SERT methylation 
levels within the CpG island by means of a partial principal 
component analysis, which overall revealed no significant results 
(Supplementary Information 3). 

DISCUSSION 

In the light of overall conflicting findings on whether the 5-HTTLPR 
S allele conveys disease vulnerability upon environmental 
adversity, 7,8 this study aimed to explore molecular modifications, 
which possibly mediate GxE interaction. Here, we report that both 
the S allele and prenatal/early adversity associate with decreased 
peripheral SERT mRNA levels in an additive manner and, 
remarkably, account for a comparable amount of variance in 
SERT expression. These effects appeared to be largely indepen- 
dent of methylation profiles within the SLC6A4 promoter- 
associated CpG island. 

Our finding of lower SERT mRNA levels in S allele carriers closely 
parallels previous data from in vitro research 2,1 9-21 but stands at 
variance with several brain imaging studies reporting mixed 
results regarding allele-specific SERT availability 24-26 Besides 
inconsistencies related to methodological aspects (see Willeit 
etal. 62 ), SERT mRNA may simply reflect a more proximate measure 
of transcriptional activity than SERT availability. Indeed, research 
across different species has shown SERT mRNA to be subjected to 
complex post-transcriptional regulation, such as translational 
repression by miRNA. 63-65 In addition, our findings implicate that 
differential exposure to early adversity may either overshadow or 
pronounce allele-specific effects on in vivo SERT mRNA expression 
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Table 2. Main effects of 5-HTTLPR genotype and stress/trauma-related variables on SERT mRNA expression profiles and mean SERT methylation 
levels (mean ± s.d.) 








Group factor: 5-HTTLPR genotype 








Additive model (tri-allelic) 




5 allele dominant model (tri-allelic) Total 




C /C in) 

5/5 (n = 18) 


SL (n = 68) 


LL (n=47) 


P 


55/5L (n = 86) 


LL(n=47) P 


SERT mRNA expression 
(relative quantification) 


0.051 ±0.03 
(0.058 ±0.03) 


0.059 ±0.05 
(0.058 ±0.05) 


0.079 ±0.06 
(0.085 ±0.07) 


0.07 
(0.04*) 


0.057 ±0.05 
(0.058 ±0.05) 


0.079 ±0.06 0.03* 0.06 ±0.05 
(0.084 ±0.07) (0.01**) 


iviean otm metnyiation [yo) 


5.1 9 ±0.62 
(5.1 6 ±0.74) 


5.08 ±0.80 
(5.04 ±0.77) 


5.04 ±0.77 
(5.1 4 ±0.80) 


0.80 
(0.75) 


5.11 ±0.77 
5.07 ±0.77 


5.05 ±0.77 0.50 5.09 ±0.77 
(5.1 4 ±0.80) (0.73) 








Group factor: trauma/stress 






Prenatal trauma/stress 




Early trauma/stress 


Recent trauma/stress 




Yes (n = 39) 


No (n = 46) 


P Yes (n = 


17) No (n 


= 116) P 


Yes(n=34) No (n = 99) P 


SERT mRNA expression 
(relative quantification) 


0.049 ±0.04 


0.073 ±0.05 


0.02* 0.041 ±0.01 0.069 ±0.01 0.04* 


0.061 ± 0.05 0.075 ± 0.06 0.1 9 


Mean SERT methylation (%) 


5.05 ±0.81 


5.05 ±0.67 


0.68 5.07 ±0.77 5.09 ±0.78 0.86 


5.1 2 ±0.80 4.99 ±0.65 0.39 


Abbreviation: SERT, serotonin transporter. Asterisks indicate significant differences between groups (*P<0.05; ** 


P<0.01). 



(Figure 2) and may thus have contributed to previously observed 
discrepancies. 22,23 

To date, the link between environmental adversity and 
persistent changes in SERT mRNA profiles has almost exclusively 
been addressed in non-human research. Specifically, maternal 
separation has been associated with lower raphe SERT mRNA 
levels in rodents 28,29 and decreased SERT availability in non- 
human primates. 66 Likewise, Rhesus macaques exposed to 
maternal aggression were found to have reduced SERT mRNA 
levels in peripheral blood cells, indicating that stress-induced 
changes of SERT expression are not limited to the brain. 31 Our 
study complements and extends previous animal research by 
demonstrating long-term signatures of early adversity using easily 
accessible markers of human SERT expression. Here, we provide 
first evidence for reduced SERT mRNA levels in individuals 
exposed to maternal prenatal stress or child maltreatment. This 
observation is strengthened by an inverse correlation between 
SERT mRNA levels and the magnitude of prenatal/early adversity 
within the overall sample. Strikingly, we observed no such effect 
for adult stressors/trauma, hinting towards a sensitive period in 
early development 67 

Epigenetic modifications are considered a promising pathway 
mediating sustained changes of gene expression in response to 
early adversity. 32 Methylation profiles within a 799-bp CpG island 
in SLC6A4 have recently been associated with SERT transcription 
and are further responsive to environmental influences 68 Several 
in vivo and in vitro studies found site-specific SERT methylation to 
promote gene silencing in transformed lymphoblast cell lines, 23,33 
peripheral blood mononucleated cells 35 and buccal cells, 34 in 
some studies depending on 5-HTTLPR. 40 Our results tentatively 
concur with previous observations by demonstrating negative 
correlations of SERT methylation and mRNA expression for 10 out 
of 83 CpG sites investigated. Although this association remained 
significant only for CpG9 after Bonferroni correction, the finding of 
12% of CpG sites being functionally relevant on an uncorrected 
level is unlikely to result from chance 10|a = 0.05, n = 83) 
< 1%). Despite its putative role in transcriptional regulation, we 



found very limited evidence for SERT methylation mediating 
associations between environmental adversity and SERT mRNA 
expression. Regarding childhood trauma, our study conflicts with 
previous reports linking site-specific SERT methylation to a history 
of sexual abuse, 23,36,37 childhood trauma 39 and bullying 
victimization. 38 However, it is of note that no specific CpG site 
has yet been consistently associated with early adversity across 
previous studies. Variable findings may result from diversity of 
used cell populations, examined subregions, methylation detec- 
tion methods (for example, quantitative mass spectroscopy, 
Sequenom EpiTYPER platform, pyrosequencing), type of trauma 
and lack of correction for multiple testing. 

Regarding prenatal stress, initial evidence has suggested a 
positive association between maternal depressed mood during 
the second trimester and the mean methylation within a 
subregion of the SERT CpG island (10 CpG sites). 69 Whereas we 
observed no effect on the mean SERT methylation, individuals 
exposed to maternal prenatal stress were characterized by 
increased methylation levels at four CpG sites (Bonferroni- 
corrected at CpG30 and CpG9, at trend level). Although it is 
tempting to suggest CpG9 methylation as a putative mediator of 
lower SERT mRNA levels in response to prenatal adversity, caution 
is advised when interpreting this finding. A closer inspection of 
CpG9 (Figure 3b), but also of some relevant sites identified by 
previous studies, revealed that absolute methylation appears to be 
marginal, falling below the detection limit for the majority of 
individuals, and is thus unlikely to constitute a valid candidate. 
Against this background, candidate SERT CpG sites should be 
carefully selected in consideration of previously observed absolute 
methylation levels in future studies. Elucidating the precise 
mechanism of CpG site-specific methylation mediating gene 
silencing, such as altered transcription factor binding, might 
further advance this selection process. The lack of robust findings 
does not rule out the possibility that methylation patterns outside 
of the investigated region or epigenetic modifications other than 
methylation may have mediated effects of early adversity on SERT 
expression. Interestingly, genome-wide analyses have suggested 
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Figure 2. SERT mRNA expression levels (mean ± s.e.m.) calculated with the 2 2AACt method as a function of 5-HTTLPR genotype and (a) 
prenatal stress or (b) early trauma. Asterisks indicate significant differences between groups (*P<0.05). 
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4741 gtcccaggaagaaagagagagcagctt tcg' ggatggggacg^atggggaggtgtccg^ggt 

4801 caagagaaagcg 6 gcacg 7 agcagacccctgtgtgccg H tcctgtggg i gjr2cg 10 gggcg 11 gcaggg 

4861 gaggcg 12 cacacctgctcctttgtgcagcctcccccctcccg 13 caaagtaaagagcaggaa 

4921 agtcaggattcctcg 14 ctcg 15 gccctgccctgccg 16 gctgctccg 17 cg 18 ctccg 19 ctcctccctgc 

4981 g 20 agcg 21 tgtgtgtgtgtcg"ggggtccctcccctcctggctctggggtcg 23 ggcg 24 c2^cacccc 

5041 g 26 cccg^tag cg 28 cg^ gcccctccctggcg 30 agcg"caaccccatccagcg^ggagcg"cg 34 gagccg 35 

5101 cg 36 gccg 37 cg 38 gggaagcattaagtttattcg 39 cctcaaagtgacg 40 caaaaattcttcaagagc 

5161 tctttggcg 41 gcg 42 gctatctagagatcagaccatgtgagggcccg 43 cg 44 ggtacaaatacg 45 gc 
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5341 tctccg 62 cctcg 63 aacccacg 64 tttctttccagaccttcttccccg 6S cctcg 66 gggagggggata 

5401 gaaccg 67 ctgcg 68 ccccaccg 69 ccctgcg 70 aggagg cg 71 aggaggtgcatgcg 7 ^ccccagcq 73 gtgg 

5461 gcg 74 ccg 75 gatcctgcccctgcg 76 ccctccacg 77 ctcagcaagagccagagctgaagctgaccg 78 

5521 gccagagtgggagacg ,9 aggaacg 80 tggagtgctcg al aagtgggcg 82 gg cg 83 tagggggctcctt 




CpG Site 

Figure 3. (a) The sequence of the SERT promoter-associated CpG island as previously described by Philibert et al. 33 (GenBank accession 
number: NG_01 1747). CpG sites analyzed by means of bisulfite pyrosequencing are numbered and base pair positions according to the NCBI 
genome browser are depicted on the left hand side. CpG sites associated with reduced SERT expression are marked by single underlining 
(uncorrected significance) and double underlining (Bonferroni-corrected significance). CpG sites with increased methylation levels in 
individuals with prenatal stress exposure are marked by light shading (uncorrected significance) and dark shading (Bonferroni-corrected 
significance), (b) Boxplots displaying DNA methylation levels across the 83 CpG sites. The box covers the methylation data of each CpG site 
between 25th bis 75th quantile (median ± one interquartile range), the whiskers represent the range of values that fall within 1.5-fold the 
interquartile range. The horizontal line reflects the methylation detection limit for bisulfite pyrosequencing. 



the possibility that methylation patterns associated with gene 
expression are more likely to be located outside CpG islands. 70 
Supporting this notion, a recent study reports that methylation in 
the shore of the SERT CpG island upstream from exon 1A predicts 
SERT expression, thus identifying this region as a potential target 
for future studies. 23 

Several limitations of the present study should be acknowl- 
edged. First, our findings rely on retrospective self-report 
measures of environmental adversity, which could be subject to 
bias. However, the finding of reduced SERT mRNA expression 
being a correlate of early adversity appears to be consistent across 
different sources of self-report. Furthermore, larger studies are 
needed to evaluate differential effects of specific types of 
stressors/trauma on SERT expression levels. Second, findings of 
the present study that were obtained in a homogeneous sample 
of healthy young Caucasian individuals may not generalize to 
other ethnic groups. Third, it remains to be elucidated whether 
findings obtained with peripheral markers of SERT expression and 
methylation profiles generalize to brain tissue. Despite being 



tissue-specific, post-mortem studies have revealed substantial 
correlations across peripheral and neural cells for both gene 
expression 71 and DNA methylation patterns. 72 Regarding SERT in 
particular, a recent brain imaging study 73 indicates that peripheral 
SERT methylation is indeed an informative marker for in vivo 
serotonin synthesis. Even more importantly, environmentally 
induced modifications in gene expression and methylation 
patterns appear to be system-wide 44,45 supporting the usefulness 
of peripheral markers for the present study. Lastly, the choice of 
analyzing whole blood was motivated by the fact that it is readily 
available, allows to stabilize RNA at the time of blood draw and, 
importantly, does not require transformation known to modify 
methylation profiles. 74 Despite these advantages, the heteroge- 
neous mixture of cell types may constitute a potential confound, 
although expression levels of very few mRNA transcripts were 
found to depend on cell composition 47 

In conclusion, our findings raise the assumption of lower SERT 
mRNA expression being a correlate of both the 5-HTTLPR S allele 
and prenatal/early adversity. Strikingly, individuals carrying both 
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the genetic and one of the environmental risk factors were found 
to have 32.8% (maternal prenatal stress) and 56.3% (child 
maltreatment) lower SERT mRNA levels compared with those 
without any risk factor. Sustained alterations in SERT expression 
profiles may in turn relate to changes in serotonergic functioning 
and thus constitute a possible molecular pathway mediating 
disease vulnerability as a function of GxE. In accordance with 
genome-wide transcriptomics analyses, 46,47 preliminary evidence 
from our pilot study suggests that interindividual variability in 
SERT expression is substantially higher compared with intraindi- 
vidual variation, thus highlighting the potential usefulness of this 
specific transcript as a biomarker. Although not without 
inconsistencies, 75-77 numerous studies indeed suggest that lower 
SERT expression in the brain and periphery associates with major 
depression, 14-16 treatment efficacy 1 4,78,79 and increased stress 
sensitivity. 17,18 Since the mean age of onset for depression 80 is 
slightly higher compared with the mean age of our sample, 
altered SERT expression levels could be discussed in terms of a 
premorbid risk factor for the development of stress-related 
psychopathology later in life. Future evaluation of such putative 
biomarkers of disease vulnerability in longitudinal designs may 
not only shed light on the pathogenesis of stress-related 
psychopathology but may also bear important implications for 
predictive testing. In line with this assumption, peripheral blood 
gene expression profiles measured in the third trimester of 
pregnancy have been found to reliably predict subsequent 
postpartum depression within a high-risk cohort. 43 Moreover, 
peripheral markers of gene expression may help to identify unique 
disease-related transcriptional signatures associated with specific 
genetic and environmental risk factors as shown for post- 
traumatic stress disorder 42 
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